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SUMMARY 


A solution to the problem of finding the energy -optimal trajectory for 
terminal rendezvous with a satellite station in circular orbit by using the 
classical calculus of variations is described. A terminal -stage engine with 
constant thrust and fixed burning time is specified, and both the optimal 
initial conditions and thrust orientation time history are determined by 
appealing to the criterion that the energy gain be maximized. The problem 
is formulated in three dimensions, with the constraint that the terminal 
stage must turn through a prescribed initial angle between the target orbital 
plane and the rendezvous -vehicle velocity vector. 

Numerical solutions to the two-point boundary-value problem were gener- 
ated by using an adjoint iteration technique. Near-earth orbital rendezvous 
trajectories were calculated for orbital altitudes of 100 , 200 , and 500 nau- 
tical miles, initial thrust-to-mass ratios of 0 . 25 g and lg, and burning 
times of 100 and 200 seconds. The thrust angles for these solutions are 
found to lie very close to the local horizontal in elevation and to maintain 
a nearly constant bearing in azimuth with respect to the orbital plane of 
the station. 

Quasi -optimal trajectories, for which the thrust vector is constrained 
to the local horizontal plane and to a constant azimuth angle with respect 
to the station plane, were computed. These trajectories compare very 
favorably with the optimal case. A degradation in effective characteristic 
velocity of only 0. 063 percent is shown in the worst case (high thrust, long 
burning time, low altitude, and maximum plane-change angle). Gravity -turn 
trajectories, which were computed for the coplanar case only, proved to be 
slightly less efficient than the qua si -optimal trajectories in each case. 

The maximum coplanar effective characteristic velocity degradation is 0. 0^0 per- 
cent for the quasi-optimal mode, as compared with a 0. 046-percent degradation 
for the gravity-turn case. 



IMRODUCT ION 


Studies of trajectories for satellite rendezvous and of the associated 
guidance problem have appeared frequently in the literature of the past few 
years. A summary of a large sector of this work is presented in reference 1, 
which is particularly valuable for its extensive bibliography. Typical 
studies directed toward satellite rendezvous trajectories are found in 
reference 2 where coasting trajectories for various incremental velocity 
changes at rendezvous are examined. 

A logical extension of previous work is consideration of the finite 
thrust problem, and, in particular, the generation of optimal thrusting tra- 
jectories. Reference 3 uses bounded thrust and a coast to achieve a point- 
to-point transfer in a specified time for examining optimization of the 
coplanar rendezvous problem. This approach treats the linearized equations 
of motion in an attempt to develop a guidance scheme amenable to onboard 
computation. 

The current paper contains the results of an investigation of terminal- 
stage guidance in three dimensions for rendezvous with a satellite station 
in circular orbit. The particular task was one of determining the initial 
conditions and steering for a constant -thrust vehicle with fixed burning time 
such that the total energy gained during the terminal stage is maximized, 
and a desired change of the orbital plane is achieved. 

Classical techniques of the calculus of variations were applied to 
determine the optimal steering function, and an adjoint iteration technique 
was employed to yield numerical solutions. Rendezvous with a near-earth 
satellite was chosen as a significant problem to demonstrate the method. For 
this problem, parametric effects were investigated by varying the ratio of 
thrust to initial mass, the burning time, the station orbital altitude, and 
the required angular change of the orbital plane. The values of these 
parameters were: initial thrust -to -mass ratios of 0. 25g and lg; burning 
times of 100 and 200 seconds; station altitudes of 100, 200, and 500 nautical 
miles; and plane changes up to the limiting capability for each combination 
of thrust and burning time. 

The use of analytical guidance schemes which provide reasonable approxi- 
mations to the optimal solutions is generally more practical than attempting 
to mechanize the exact optimal steering scheme. The basis of such schemes 
is a simplified steering law to specify the thrust direction. Two simple 
laws were investigated. The first, called the quasi-optimal case, was derived 
from inspection of the optimal solutions. This scheme consisted of con- 
straining the thrust vector to the local horizontal and holding a constant 
bearing angle with respect to the satellite orbital plane. The second scheme 
was the gravity turn, where the thrust vector is always alined with the 
velocity vector. These laws were evaluated by comparing their energy gain 
with that of the optimal solution. 



The author wishes to acknowledge the assistance of Mrs. Beverly P. 
Latimer and Miss Harriette A. Seals of the Analysis and Computation Division, 
NASA Langley Research Center, who programed and checked the digital computer 
mechanization of this problem. 


SYMBOLS 


A 

B 

* 

c 

E 

F 

f 

g 

h 

s 

I 

K 

m 

m(0) 

n 

P 

r 

e 


r 

s 

T 

t 

u 


coefficient matrix for adjoint differential equation set 

transformation matrix relating initial value perturbations 
to terminal errors 

effective rocket exhaust velocity, ft /sec 

2 2 

total energy per unit mass, ft /sec 
modified function to he optimized 
arbitrary function 
acceleration of gravity, ft /sec 
satellite station altitude, nautical miles 

integral 

constant 

mass of vehicle, slugs 
initial mass of vehicle, slugs 
iteration cycle number 
adjoint variable 
radius of earth, ft 

rendezvous vehicle radius from center of earth, ft 

station radius from center of earth, ft 

rocket thrust, lb 
time, sec 

state variable in expanded system 


V 


total velocity of vehicle, ft/sec 


characteristic velocity, ft/sec 
station velocity, ft/sec 

X,Y,Z coordinates of rotating space station centered axes (fig. l), ft 

a thrust elevation angle with respect to local horizontal, radians 

3 angle between rendezvous vehicle velocity vector and station 

velocity, in the station horizontal plane, deg 

2 3/2 

X gravitational constant, gr^ , ft /sec 

2 2 

AE change in energy during rendezvous stage, ft /sec 

AV effective characteristic velocity, including potential energy 

contribution, ft/sec 

AV differential effective characteristic velocity, ft/sec 

Su ' perturbation on state variable u 

0 thrust elevation angle with respect to satellite reference axes, 

radians 

A Lagrange multiplier, ft 

t burning time, sec 

cp constraint parameter 

^ thrust azimuth angle, radians 

rate of station rotation, radians/sec 

Subscripts: 

i,j,k state and adjoint variable indices 
max maximum 

opt optimal 

quasi -opt quasi -optimal 

1,2,3 Lagrange multiplier indices; constraint indices 
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Operators: 

( ) differentiation with respect to time 

( ) second differentiation with respect to time 

( ) boundary value specified by problem constraints 

$( ) variational operator 


ANALYSIS 


Formulation of the Problem 

The mathematical model employed for this investigation of satellite 
rendezvous was a mass-particle with three degrees of freedom moving in a 
rotating cartesian axis system oriented in a satellite station. The station 
was assumed to be in a circular orbit about a uniformly spherical earth. 
These axes and the associated notation are illustrated in figure 1. 



Figure 1.- Axis system and angle definitions. 
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The equations of motion for this model are 


where 


and 




(i) 


A rendezvous vehicle with constant thrust directed along one vehicle axis 
and a constant burning rate was assumed. The angles 0 and r are the 
control variables and may be chosen arbitrarily. Reference 2 contains a 
derivation of these equations of motion. 

The quantity selected for optimization was the change in total energy 
over the terminal stage. The total energy is expressed as 



( 2 ) 


This criterion was used rather than the familiar minimization of expended 
mass in order that the burning time of the fixed-thrust engine could be 
specified and all the initial conditions left free. The physical interpre- 
tation of the criterion is best illustrated by an example. Consider the 
problem of rendezvous from a parking orbit consisting of an initial boost, a 
coast period, and a terminal boost. The energy (eq. (2)) remains constant 
over any coasting orbit. Since the energy of the satellite station orbit is 
fixed, maximization of the change in energy over the terminal stage leads to 
the lowest possible energy for the intermediate coast and hence to the minimum 
energy cut-off condition for launch boost. 
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The geometrical boundary conditions for the problem were then left free 
at the initiation of terminal-stage burning except for the specification of 
the angle between the orbital plane of the station and the coasting ascent 
plane of the rendezvous vehicle. It was required that the origin of the 
station-centered coordinates be reached with no relative velocity remaining 
at the specified time of burnout. 

To force satisfaction of the geometrical end constraints, the problem 
was inverted and run backward in time starting from the origin. The boundary 
condition on plane change can then be expressed as 


tan $( 0 ) =5 


X + n s( r s - Z )Lo 


( 3 ) 


For the backward problem, the optimization task becomes one of mini- 
mizing the initial total energy. The integral corresponding to the energy, 
together with the equations of motion (eq. (l)) and the boundary condition 
(eq. ( 3 )) adjoined as constraints, was varied to determine the Euler -Lagrange 
equations. Appendix A contains a detailed description of this process by 
which the equations for the control variables, the Lagrange multipliers, and 
their associated boundary values are determined. 

Trajectories which satisfy the Euler -Lagrange equations and match the 
boundary conditions correspond to relative minimums for the optimization 
problem, neglecting questions of existence and uniqueness of the numerical 
solution. In order to determine the global minimum solution within the 
stated constraints, all possible relative minimums must be rigorously in- 
vestigated. Similarly, to determine the validity of the numerical solution, 
the effect of perturbations on its stability must be determined. For the 
current problem, the axis system employed affords a well-conditioned set of 
equations of motion. Judicious search for alternate relative minimums, 
together with engineering judgment, was used to determine that the results 
represented global minimums. 

Iteration by the Adjoint Method 

Once the differential equations governing the variational problem have 
been written and the requisite number of boundary conditions has been speci- 
fied, the problem is defined in the mathematical sense. Unfortunately, 
equations (l) are nonlinear and must be integrated numerically. To make 
this problem computable, the two-point boundary-value problem must be con- 
verted to an initial -value problem. A standard approach to this process 
is to estimate initial values for the unknown boundary conditions, integrate 
over the interval, and then observe the errors in the known end values and 
operate on them to improve the initial estimates. The logical process by 
which the end errors are converted to initial-value corrections should force 
convergence to the desired end values within a reasonab.l e number of itera- 
tions. 
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For this study, an iteration process described in reference 4 was 
employed. This technique has also been employed recently in reference 5 
for the solution of variational problems which are by nature two-point 
boundary -value problems. The iteration process, together with the resulting 
equations for application to this problem, is derived in appendix B. 


RESULTS AND DISCUSSION 


Range of Parameters 

Numerical results were generated for the problem of rendezvous with a 
near-earth satellite in circular orbit. The rendezvous vehicle was assumed 
to have constant thrust directed along one vehicle axis and a specified 
thrust-time duration. A constant mass rate was used which corresponded to 
an effective exhaust velocity of 10,000 feet per second (specific impulse 
of 310 sec). 

The ranges of parameters chosen for investigation were initial thrust - 
to-mass ratio, burning -time duration, orbital altitude, and plane-change 
angle. The values used in this investigation are tabulated as follows: 


T/m(0), g 

t, sec 

3(0), deg 

0. 25 

100 

0 

.5 

1.0 

i.5 


200 

0 

1.0 

2.0 

5.0 

1.0 

100 

0 

2.0 
4. o 
6.0 


200 

0 

4.0 

8.0 
10. 0 


All cases were run for orbital altitudes of 100 , 200 , and 500 nautical miles. 
To determine the range of the plane-change angle, an upper bound in 3 was 
calculated from the following expression: 


3 


max 



(4) 


8 




where 


V 

c 


,* 



There is a different characteristic velocity for each thrust and burning - 

time combination. 


Adjoint Iteration Check 


The combined variational problem and adjoint iteration scheme were tested 
and found to be rapidity convergent even for very poor guesses of the initial 
multiplier values. Table I illustrates the convergence process for two cases. 
The conditions of case 1 are: T/m(0) = 0.25g, burning time of 100 sec, 

1 ° 

station altitude of 100 nautical miles, and a plane-change angle of 1— . The 


sum of the errors in end conditions 



drops seven orders of magnitude 


i 

in four passes. Case 2 conditions are T/m(0) = lg, burning time of 200 sec, 
and a plane-change angle of 10°. This more severe case required 12 passes, 
but the error is seen to be uniformly convergent and no dif f iculties occurred 
during iteration. Case 2 displayed the slowest convergence of all the cases 
computed. 


Optimal Trajectories 

The profiles .of the optimum trajectories, together with time histories 
of the thrust vector orientation, are presented in figures 2 to 5* The 
origins of the curves of the thrust elevation angle are displaced for clarity. 
The thrust elevation angle measured with respect to the local horizontal a 
exhibits variations with time and required angular change of plane. However, 
the maximum departure from the local horizontal encountered is less than 
0.055 radian. The thrust azimuth angle is nearly constant for all the cases 
run. 


The results from the optimal trajectories were used to evaluate the 
energy decrement resulting from use of finite burning times as compared with 
impulsive velocity changes. In order to make a direct comparison, effective 
characteristic velocities were calculated from the energy increment realized 
by using the following relation: 


AV 

c 



(5) 


The variation of effective characteristic velocity and of percentage 
velocity differential compared with an ideal single impulse with plane change 
and altitude is shown in figure 6. As would be expected, the energy change 
realized decreases with increasing change of orbital plane. This decrease . 












Lateral velocity, Radial velocity, Longitudinal velocity, Lateral displacement, Radial displacement, 








table i. - adjoint iteration 


Case 

Cycle 

^(t) 

>^(t) 

^(t) 

X 2 (t) 

A^(t) 

*j(t) 

AI(0) 

^(0) 

^ x (o) 

^,(0) 

1 

1 

0 

0 

0 

0 

0 

0 

-0. 2165x1c) 5 

O.2687XIO 2 

0. 7932X10 1 

0 


2 

770.5 

- 7.940 

7,127 

1.197 

46.32 

- .9176 

- .5436X10 1 

. 6105x10^ 

-.3094 

-.2036X10 1 


3 

678.5 

- 7. 633 

7,244 

3.328 

36.19 

- .7944 

. 9569X10" 2 

. 4607X10 1 

. 1020X10* 1 

. 3459X10 -1 


4 

675.0 

- 7.645 

7,246 

3.292 

35.84 

- .7923 

- .9460x10"^ 

-. 377JX10 -1 

3789X10“^ 

1348xio“ 2 

2 

1 

0 

0 

0 

0 

0 

0 

ii 

^>.3651x10 

0. 1328X10 4 

0. 2996x10^ 

0 


2 

4,-846 

-30. 62 

17,733 

15.14 

643.4 

-5.687 

- .7745X10 5 

. 2946xl0 4 

-. 7171X10 1 

-. 1107X10 2 


3 

429 

-23.32 

26,984 

31.18 

- 324.5 

- .8022 

- .3099X10 5 

. 2362x10 

-.5812 

-. 7151 


4 

-1,940 

-23.76 

34,8l8 

34.56 

-882.3 

I.881 

- . 1125x1 o 5 

. 1391X10 4 

-.1536X10 1 

-.1309X10 1 


5 

-3,567 

-22.97 

39,397 

37. 31 

-1,266 

3.879 

- . 3531X10 2 

. 86UX10 5 

-.6680 

-.5354 


6 

-4,516 

-22.78 

41,631 

38.52 

-1,488 

4.998 

- .9836XIO 1 

. 4158X10 5 

-.2489 

-. 1968 


7 

-4,960 

-22.75 

42,528 

38.97 

-1,590 

5.493 

- . 2 52 2X1 0 1 

. 1665XIO 5 

-.7707X10" 1 

62 3 3x10 ~ 1 


8 

-5,134 

-22.76 

42,843 

39.12 

-1,630 

5.677 

- .6174 

. 5842X10 2 

-.2l66xi0 _1 

-, I8l6xl0 -i 


9 

-5,194 

-22.77 

42,946 

39.17 

-1,643 

5-738 

- .1472 

. 1877X10 2 

-.5823xio“ 2 

-. 5068X10 -2 


10 

-5,213 

-22.77 

42,977 

39.18 

-1,648 

5.757 

- .3442X10" 1 

. 5686X10 1 

-. 1531X10' 2 

-. 1379X10" 2 


11 

-5,219 

- 22.78 

42,986 

39.19 

-1,649 

5.762 

- .7843X10" 2 

. 1653X10 1 

-.3964xl0“ 5 

-. 3679X10" 5 


12 

-5,22° 

-22.78 1 

42,988 | 

39.19 

-1,649 

5.764 

- . 1831XIO”* 2 

.4656 

-. 1005X10 ’ 5 

-. 9680x10 


is less marked at higher altitudes, primarily because a smaller addition of 
velocity normal to the station plane is required to rotate the orbital 
velocity through the desired angle (orbital velocity decreases with altitude). 
The ratio of characteristic velocities similarly decreases vith increasing 
orbital-plane change, which indicates that the cost of finite burning time 
increases vith the turning angle, even when normalized by the characteristic 
velocity requirement. 

For the cost criterion, finite burning times have a small advantage 
over impulsive velocity corrections in the coplanar case vith the initial 
point unconstrained if only a terminal impulse is allowed. This is not 
totally •unexpected, since, in general, two impulses would be required to 
achieve rendezvous from arbitrary initial conditions, and the sum of the two 
velocity increments would be less than the single increment for an equal 
change in energy. A heuristic explanation for this statement can be drawn 
from consideration of the rate of change of energy 


lb 




E = V • 


V + 



f 


(6) 



CONVERGENCE EXAMPLES 


4\ 5 (0) 

^ 5 (o) 

Error 

£*S s Ct) 

1 

e i (t) 

8\(0) 

6\(0) 

6^(0) 

6^(0) 

B?v 5 (0) 

5^(0) 

-0. 5000X10 2 

- 0. 9915 

0. 50153X10 5 

-0. 3466x10? 

- 0 . 7705 x 1 c 5 

0. 7940X10 1 

4 

- 0 . 7127 x 10 

-1- 197 

-46 . 32 

0. 9176 

. 1529X10 1 

- .^955X10" 1 

. 37632X10 5 

- .3459X10? 

. 9198X10 2 

3064 

- . Il66xl0 5 

-2. 130 

10. 13 

1232 

. 1380X10" 1 

- .4592X10" 2 

. 2 12 3X1 0 2 

- . 3458 x 10 ? 

. 3529X10 1 

. 1127X10 -1 

- .2583X10 1 

. 03538 

.3502 

-. 002109 

. 6820 x 10 ~ 2 

. 1797X10~ 2 

. 1472X10 -2 

- .3458x10? 







-0. 9944X10 5 

- 0 . 9325 XIO 1 

0 . 1608 x 10 ® 

-0. 442CX10 9 

-0. 4846x10^ 

0. 3061X10 2 

- 0 . 17733 x 10 5 

-0. 1514X10 2 

- 0.6434xi0 5 

0 . 5687 XIO 1 

. 2289 XIO 5 

- . 23^9X10 1 

-9353X10 T 

- .4l42xi0? 

. 4416x10^ 

-.7291X10 1 

4 

- .9251X10 

- .1603X10 2 

. 9679 XIO 3 

4885X10 1 

1572X10 2 

- . 2926 XIO 1 

. 5674X10 7 

- . 401CX10? 

. 2369 x 10 ^ 

.4321 

- .7833X10 4 

- .3386X10 1 

. 5577X10 3 

-. 2683 XIO 1 

. 3366 XIO 2 

- . 1435 x 1 0 1 

. 1949 x 10 ^ 

- .3935x10? 

. 1627X10 4 

7897 

- .4580X10^ 

- .2750X10 1 

. 3833 XIO 5 

1998 XIO 1 

. 2l48xi0 2 

- .8573 

. 7432X10 6 

- .3900X10 9 

. 9492X10 5 

1911 

4 

- .2233X10 

- .1203X10 1 

. 2224xl0 5 

-. 1119 XIO 1 

. 1235X10 2 

- .3924 

. 1731X10 6 

- . 3888 x 10 ? 

. 4441X10 5 

-. 2388 X 10' 1 

- . 8971 x 10 ^ 

- .4562 

. 1026X10 3 

-. 4953 

. 5444X10 1 

- .1510 

. 2777X10 5 

- . 3885 XIO 9 

. 1736X10 5 

.1009X10 -1 

- .3160X10 5 

- .1517 

. 3947X10 2 

1838 

. 2Q27X10 1 

- . 5 l 1 < 8 xio~ 1 

. 3417X10 4 

- .3884x10? 

. 599QX10 2 

. 8058 x 1 o' 2 

- . 1020X10 5 

- .4653X10 -1 

. 1344X10 2 

6 O 86 XIO ” 1 

.6751 

- . 162 IXIO -1 

. 3528X10 5 

- .3884x10? 

. 1904X10 2 

. 3532X10" 2 

- .3099X10 2 

- .1359X10” 1 

. 4231X10 1 

-. 1877X10 -1 

.2093 

- .4844xlO~ 2 

. 3238 XIO 2 

- .3884x10? 

. 5723X10 1 

. 1262 * 1 0 -2 

- . 903 CX 10 1 

- . 384lxi0“ 2 

. 1263 XIO 1 

-. 5524X10" 2 

. 6187X10 " 1 

- . 1393X10 " 2 

. 2736X10 1 

- .3884x10? 

. 1654X10 1 

. 4071X10 -5 

- .2549X10 1 

- . 1058X10" 2 

.3633 

-. 1571X10" 2 

. 1765X10 -1 

- .3895X10 " 5 

.2171 

- .3884x10? 








It is clear that the incremental change in energy is maximized if the vehicle 
is accelerated when the velocity is as high as possible (and with acceleration 
principally along the velocity vector). Hence for this problem, it is advan- 
tageous to thrust before the vehicle reaches apogee (minimum velocity) as long 
as no plane change is required. In the latter case, the advantage of turning 
a smaller velocity vector negates the benefits of accelerating at high veloc- 
ity. 


From inspection of characteristics of the optimal trajectories, a 
quasi-optimal flight was defined for which the thrust vector was constrained 
to the local horizontal plane and held at a constant bearing angle with 
respect to the satellite orbital plane. It was necessary to iterate to 
determine the bearing angle necessary to achieve a prescribed plane change. 

Figure 7 is a typical plot of the losses in effective characteristic 
velocity in the quasi-optimal case. The parameter plotted is differential 
effective characteristic velocity, defined as 


AV = AV 

c, quasi -opt 


AV 


c,opt 


(7) 
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Figure 6 .- 


Effective characteristic velocity variation with orbital altitude 

and inclination. 


(a) 


T 

m(0) 


0.25g; t= 100 sec. (b) 


0.25g; t= 200 sec. 


A negative AV indicates that the energy gain over the quasi-optimal tra- 
jectory is less than that over the optimal trajectory. Note that the losses 
decrease and then increase with increases in orbital -plane change. This 
variation can he correlated with the time histories of optimal thrust eleva- 
tion angle in figures 2 to 5 where it can he seen that this angle most closely 
approaches zero for medium orbital -plane changes. The comparison between 
optimal and quasi-optimal paths is very favorable, which indicates that 
guidance modes based on the quasi-optimal path should be highly efficient. 

The maximum penalty for using the quasi-optimal guidance mode occurs for the 
maximum plane change and is only 0 . O 63 percent for the worst case encountered 
(T/m(0) a lg, t = 200 sec, h g s 100 nautical miles). 

A second simplified flight path, the gravity -turn trajectory, was tested 
for the coplanar cases. 

In the gravity turn, the thrust is always collinear with the velocity 
vector which lies close to the local horizontal for rendezvous with a 

16 




Plane change angle, P, deg Plane change angle, P, deg 


Figure 6. - Concluded. 

(°) jujo) = lg j T= 100 sec * m^o) = lg > T= 200 seCl 


satellite* in circular orbit when the velocity gains over the terminal stage 
are of the order considered here. The resultant effective characteristic 
velocities vere calculated by using equation (5) and are listed in table II, 
together with comparable figures for the optimal and quasi -optimal cases. 

The quasi-optimal mode proved to be slightly more efficient than the 
gravity -turn mode. Effective characteristic velocity showed a maximum 
degradation of 0. 040 percent for the case where T/m(0) = Ig, t = 200 sec, 
and h g = 100 nautical miles, as compared with a 0. 046-percent degradation 

for the gravity turn. 
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Figure 7 *- Variation of effective characteristic velocity differential 
•with orbital inclination angle and orbital altitude for q_uasi-optimal 


T 

trajectories, - lg; t 


200 sec. 


TABLE II. - EFFECTIVE CHARACTERISTIC VELOCITY FOR COPLANAR CASE 


h . nautical miles 
s 

T 

JtfO? g 

t, sec 

AV . , 

c , opt 7 

AV, ft /sec 

ft/sec 

Quasi -optimal 

Gravity turn 

100 

0. 25 

100 

77*+. 26872 

-0. 00201 

-0. 00202 



200 

1,495.2660 

- .0591 

- .0592 


1.0 

100 

2,791.8851 

- .1049 

- .1101 



200 

4,977.6285 

-1.9765. 

-2.2859 

200 

0. 25 

100 

774. 26198 

-o. 00162 

-0. 00163 



200 

1,495.1625 

- ,0468 

- .0486 


1.0 

100 

2,791.7875 

- .0849 

- . 0898 



200 

4,976.2501 

-1. 9188 

-2. 1217 

500 

0.25 

100 

774. 25525 

-0. 00127 

-0. 00130 



200 

1,495.0751 

- .0575 

- .0389 


1.0 

100 

2,791.7053 

- .0682 

- .0724 



200 

4,975.0704 

-1. 6754 

-1. 8387 


18 



CONCLUDING REMARKS 


Energy-optimal trajectories for terminal-stage rendezvous with a 
satellite station in circular orhit which used a constant -thrust rocket with 
a fixed burning time have been derived by using the calculus of variations. 

The problem is three dimensional, including the constraint that the terminal 
stage is to turn through a given initial plane-change angle. Numerical 
solutions to the resulting two-point boundary-value problem have been gener- 
ated by using an adjoint iteration technique. 

Near-earth orbital rendezvous was examined for orbital altitudes of 
100 , 200 , and 500 nautical miles, initial thrust -to -mass ratios of 0 . 25 g 
and lg, and burning times of 100 and 200 seconds. The resulting prescribed 
thrust angles lie very close to the local horizontal in elevation, and are 
nearly constant with respect to the target orbital plane in azimuth. 

Quasi -optimal trajectories were computed by holding the thrust vector 
in the local horizontal at a fixed azimuth angle. These trajectories, which 
show a maximum degradation in effective characteristic velocity of 0. 065 per- 
cent, compare very favorably with the optimal trajectories. Gravity turns 
were also run for the coplanar case. 

The quasi-optimal path proved to be the more efficient of the simplified 
trajectories. A degradation in effective characteristic velocity of 
0 . 040 percent was shown for the worst coplanar case as compared with 0.046-per- 
cent degradation for the gravity turn. 


Manned Spacecraft Center, 
National Aeronautics 
Houston, Texas, 


and Space 
March 26 , 


Administration, 

1963. 
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APPENDIX A 


OPTIMIZATION OF ENERGY CHANGE FOR NONLINEAR EQUATIONS OF MOTION 


The description of the problem of rendezvous of a mass -particle with a 
satellite in circular orbit is provided by the following equations of motion: 


X = - cos 0 cos t + 2Q Z - x/— U - Q 2 ' 
in s 3 s 


(Al) 


Y = — cos 9 sin \|r - Y * 

r/ 


(A2) 


= - — sin 0 - 2QX + fr - ZV| 
m s \ s J\ 


“■ - Q 

3 s 


(A3) 


The variables are defined with respect to the satellite in a rotating system. 
(See fig. 1. ) It is required that the change in total energy over the burning 
period be maximized. The variational notation used in reference 6 is employed. 


The energy is expressed as 


v 2 r 

E -T -J- 

f 


(A4) 


where 


V 2 = X 2 + Y 2 + z 2 + 


X(r - Z 1 ) + ZX 

- Cl 2 

X 2 + (r - Z) 2 

V s ) 

s 

Vs ) 


and the integral to be maximized is 


I « / E dt 

0 


(A5) 


where 


T 

E s - 

m 


+ Q g ^r g - Z) cos 0 cos i|r + Y cos 0 sin i|r - (z + Q^X^sin 0^ 


Three constraints, and f-^, arise from requiring satisfaction of the 

dynamic equations (Al) to (A3). A fourth constraint is introduced when it 
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is required that a specified change in the orbital plane of the vehicle be 
made. This constraint can be expressed as 


<p 4 


= Y - 


X + Q It - z) 

tan 3(0) 

s\ s y 



= 0 


(A 6) 


t~0 


The first three constraints are adjoined to the function to be optimized with 
Lagrange multipliers 


F = E + A 1 cp 1 + A 2 cp 2 + 


(A7) 


The integral to be operated on is altered accordingly 

T 

r 

i 

0 


= l + hh + A 2 9 2 + ^3 cp 3) dt 


(A8) 


and the variation taken 

61, f 

0 L 


5E 


59, Scp 2 ^ Scp 5 


S + \ + A 2 + a 3 zrj 5x + [^ + ^ r 


1e . , Scp 3Uv 


- , /bw\ • 

+ A i 5x + [\ W + \ W- + A 3 SYT Y + (^) 6Y + A 2 5Y 




59, 


'59, 






5z + A i 5z + A 2 5z + A 3 bz j 5Z + (~ + A 


8Z 


v5z 


5z 


5e , , ^ + a 


59 


+ ^ 62 + (If + \ sr + \ + *s 


69 




dt 


(A9) 
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Integrating "by parts in the usual manner 


81 = 


A, SX + A 8Y + A, 8Z + f — + A, — ^ - aIsx + (~ - Aj SY 


^3X 


<*P, 

3 Si 

T A- • 


3y 




dt 


- vV 

8Y + 

W s jj 



39. 

\,Sx 5 Sx 1 

39 




5X + 


3tp, 39 39, 

A l 3F + a 2 3T + A 3 3Y 


SE ^P- 

3z + A i 3z“ + A 2 3zT + A J Sz 


^ ♦ A, ^ - h 

\3z 


dt 


1 3Z 


ai . „ , ^2 


39. 


8Z + ^ + a 


39 T A i 3T + a 2 3T T a 3 3e 


+ A, 


150 


1 3E ^1 . ^2) S\ 

+ + A i sr a 2 wrV dt 


(A10) 


Considering the boundary conditions, all the physical variables are defined 
at t = t, since it is specified that the vehicle is to back off from a 
rendezvous. Therefore, all the Lagrange multipliers are unspecified at t = t. 
At t = 0, constraint equation (A6) is applied, but the end conditions are 
otherwise free. The variation of equation (A 6) is written 


SY = (SX - Q g SZ^ tan pTo)| t=0 


(All) 


This relation is substituted to eliminate 5Y from equation (AlO), and the 
fundamental lemma of the calculus of variations is applied to the boundary- 
value bracket. This leads to the following end conditions on the multipliers: 


A 1 (0) 

= - 

A 2 

tan p(0)| t=Q 

\(o) 

11 

3 11-3 

COS 

9 cos 

1 1=0 

\(0) 

_ T 

m 

cos 

0 sin \|r| 

t=0 

a 3 (o) 

= 0 




a 5 (o) 


Q 

s 

tan p(0) - — sin 0 
m 


(A12) 
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Note t ha t A 2 (0) is free, hut Y(°) is defined "by equation (A 6). 


Again, applying the fundamental lemma to the quantities under the integral 
and thus setting the coefficient of each variation equal to zero yields the 
Euler -Lagrange equations. 

,21 ^ 


\ = 2Q A, - 

1 s 3 


- Q 

3 s 


2\ 3i2L 


r f r f 


5 »5 


+ 2 

m 


A^ = 


- 9 cos ^ s in 9 - ^ cos \|r + f sin \|/jcos 9^ 

9 sin \|/ sin 9 


„ 5 A i 13 


3^\. £1 K 


r f r f 


5 T2 


l x .2 

5 A 3 m 


r - Z 
s 


+ [ — sin ilf — ilr cos ilr^ cos 9 
\ m 


A, = 


-20 A - 
s 1 


3rx(r s - Z) ^ 3rY(r s - z 

r 5 A 1 " r 5 
f f 


N> " 


_ - 0 
3 s 


2\ > (A13) 


3r 


( r s - z i 


T 

A + — 
3 m 


— sin 9 + (q cos \|r - Q^cos 0 
m \ s Y / 


tan 9 


( Z + XQ s - A^cos jr 

X + Q (r - - A n 

s\ s / 1 


Y - A, 


tan \|r = 


X + 0 s (r s - z) - 




These equations, together with the constraint equations (Al) to (A3), the 
six boundary values at t = 0 specified by starting at rendezvous (all relative 
velocities and displacements zero), and the six boundary values at t = t from 
equations (A6) and (All), completely specify the problem. 
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APPENDIX B 


EQUATIONS FOR ADJOINT ITERATION 


A general method for the satisfaction of dynamic two-point boundary -value 
problems is described in reference 4. This method is applied to the optimi- 
zation problem of this paper. In this appendix, the iteration technique is 
first developed in general terms, and then the specific substitutions used in 
solving this problem are stated. 

The set of three second -order differential equations (eqs. (l) to (3)), 
together with the associated Lagrange multipliers (eq. (A13)) f is considered 
as the equivalent set of 12 first-order differential equations: 





u. = f . (u‘. ) 

l 3 

(Bl) 

where 

u.(0) 

-V 0 ' 5 ' 

u 

H 

V* 

ro 



^(t) 

= ^(t) 

(i = 6, 1, . . . , 11) 


The variable 

8u. is 

defined 

as a first-order perturbation about u 



Su. 

i 


I 


be. 

Su. 

OU J 
d 


(B2) 


The equations adjoint to the set of equations (B2) are described by 



(B3) 


If the set of equations (B2) is multiplied by p^ and the set of equations 
(B3) by Su^ and the resulting relations added and summed on i, then 


i 





(B4) 


It is seen that the right-hand side vanishes because of the symmetry of the 
indices. Finally, integration of equation (B4) from t = 0 to t = t yields 
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the desired adjoint relationship 

^ P 1 (0)8u i (0) = ^ P 1 (t)5u 1 (t) (B5) 

i i 

Now, if &uu(t) is defined as error observed in the end boundary values 

Su^t) = u ± (t) - u ± (t) (B 6) 

It is seen that equation (B5) gives a means of reflecting these terminal 
errors back to determine initial value improvements. 

To determine the six required values of 5u^(0), the adjoint equations 

(B3) are integrated backward from t = t to t = 0 for six distinct sets of 
end conditions. For convenience, these end conditions were chosen to be 
identically zero except for one value for each set which corresponds to a 
known end condition. 


p i,k (T) = 1 


p ± k (T) =o (k ^ i + 5 ; i » l, 2 , • . 12 ; k - i, 2 , 



(B7) 


Once this integration is performed, the resulting values are substituted into 
equation (B5 ) 


X P i,k (0) Bu i< 0) " 6u i+k (T) 


CBS) 


and the six equations yielded are solved simultaneously for 5uu(0). 
Finally, improved initial estimates are calculated 




u^(0) * u i (0) + 5^(0) 


(1 - 7 , 8 , . . 12 ) 


(B9) 


This iteration process is cycled until the end errors fall within a desired 
tolerance. 

It must he noted that the Su^(0) corrections only satisfy conditions 

for small perturbations about the nominal solution. An additional test was 
provided in the actual mechanization of the problem to prevent the computa- 
tion from diverging if the initial errors were too large. This test required 
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satisfaction of the relation 



&U . ( T ) 
1 


L U !^ 


J n+1 


&u ± (t) 

l' 


(BIO) 


n 


Nov the substitutions required to state the explicit problem are defined 


as: 


u i = x = f i 


u 2 = x = f 2 

u = Y = f. 


= Z = f 4 


u 5 = z " f 5 

"6 " X * f 6 


b' \ f 7 


u 8 " A i " f 8 


u 9 = X 2 f 9 

u l0 = X 3 f l0 
u ll = X 3 = f ll 


U 12 \ f 12 






J 


The desired boundary values on 6u. (eg. (B2)) are 


(Bll) 


&u^(t) = 0; i = 1, 2 6 
Sug(O) = [x + fi s (r g - Z) tan 0(0) - Y 


"b 


8u^(0) = - A 2 tan (3 (0 ) - ^ 


T 

8ug(0) = — cos 9 cos \|r - A., 


m 


T 

5u^(0) = — cos 0 sin i|r - A_ 
9' m 2 

8u (0) = - A I 

10 5l t =o 


t=o 


t=0 


t=0 


t=0 


&Un(°) = \ n s tan p(0) - - sin 0 - 


t=0 


J 


(B12) 


Equation (B2) is not integrated since all of its boundary values cannot 
be specified. Equation (B3) is integrated, and since it is linear and will 
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be used with six distinct sets of starting values on p., it is most easily 
treated as a vector -matrix equation. 


— — 


— 


P i,j 

J 

II 

A. V 
i,k 

P k, j 


(B13) 


where the matrix fp-, .1 is of dimension 12 X 6 and [a. | is a 12 X 12 

L k >jJ r _ L l ^ k J 

square matrix. In the matrix p L the j columns are associated with the 
six sets of end conditions. The matrix [p(t)] is defined as 


P, . ( T ) = 0 


P-j . (t ) = 0 

x > o 


p i,j (T) - 1 


(i ± 0 + 5) 

(i - 1, . . 5, 12) ) 

(i = j + 5 = 6, 7, . . .,11) 


(Bl4) 


J 


Integrate from t = 0 to t = t together with the set At 

t = t, the following relation is solved to find the required initial improve- 
ments: 

r _ r i 

(B15) 




where ( 6u t> is a six vector and [B] is a 6x6 matrix defined from the 
elements of £pj as 

M' 1 - [pj,i] < bi6 > 

where 

j = 1, 2, . . . , 6 
i = 7, 8, . . 12 

This process is equivalent to the solution of equation (B8). 

Lv] are listed below, where missing elements are 

2 


The elements of 
understood to be zero. 


\z = m Sin 0 COSt l + n' % 2 ) - 


(B17) 


I I 
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a T . . . , de 5rT£ 

\,h - - sm 0 sm * s - -y 

r f 


(B 18 ) 


T 


as . 3rx(r s - z) 


a i,6=^ cos 9 1 + — 1 


(B19) 


A, 



5X 


XA 


2 - K - z ) a 3 


5T 


Q g - 9 cos tjcos 9 + cos i|f + \|r sin i|Asin 0 


a§ 

dX 


as i . Q Svjf . Q 

- cos \|f sin 0 - ^ sm cos 0 


(B20) 


A. 


1 , 1 ° 


= _ 2r_ 

„ 5 


YAij 1 - 2 ] - XA 2 |1 2 I 


s!\ ^( r s - z ) 
2 " 2 
r f / r f 


A, 


{[• 

a© . 


T fp 

+ ~ < |^9 sin \|f cos 9 - ( “ sin - i|r cos \J/] sin 0 


a \|/ 


+ ^ sin \)r sin 0 - ^ cos i | i cos 0 


(B21) 


A. - - K {(r s - z)X,(l - 2 ^ - Z ) A, 


1,12 5 

r 

f 


11 2 
f , 



~ cos 0 - cos \|r - 9^ sin 0 


a© 

iax 


V - - 1 


(B22 ) 
(B23) 
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\ = — I 

2,2 in 1 


, T . 

2,4 mV 


ao J 

0 sin i | r 

h\f 

— + cos 

Am 

« 

ax 


ax 

a© 

0 COS 1 

a^r 

— - cos 


ax 


ax 


A 0 ^ = — cos 0 - 2Q 

2, 6 m 5x 


2,8 m 


- 0 cos ty^cos 0 + ^ cos i|f + \|r sin il^sin ©J - — 


0 sin 0 sin ^ + ( ~ sin ^ t cos \|/)cos 0 


JdX 


S0 , • Q ^ , D 

cos i|r sin 9 - — 7 sm i|r cos 0 

SX ax 


A 


T 


'2,10 m \L 


* \ ” ^0 
9 sin cos 0 - ( — sin \|/ - cos ty) sin 0 — 

j 
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0 cos \1 r sin 0+1— cos + \|r sin \|f)cos 0 

m 1 1 ax 


Sx 
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e]^ 

_ a> 


S0 




h sin \|f sin 0 - — f- cos i|r cos 0 

ax ax J 


a = _2 

12,12 m 


— cos 0 - ( D cos \|r - 0)sin oj-^r 

m ' s 1 Jax 


n sin y - cos 0 -X - — cos 9 

s ax ax 7 


(B24) 


(B25) 


(B26) 


(B27) 


(B28) 


(B29) 


29 



*3,8" 
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YA. 
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10 


2L- 
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■ 5 
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2i 


^ T 30 0 
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T / JX0 ^ , S\lr 
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SY 


d^X 
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(950) 

(951) 
(B32) 

(B33) 

(B34) 

(B35) 

(B36) 

( 937 ) 
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A, | = — (sin 0 sin ll r — - cos 0 cos f 

4 > 4 m V S Y ar/ 


n T _ 30 

A, ^ = — cos 0 — 

b ,6 m 


\8 


'"I 


(- 


- 0 cos \|r] cos 0 + ( — cos f + t sin t| sin 0 
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